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Abstract 

Until now the representation (i.e. plotting) of curve in Parallel 
Coordinates is constructed from the point ^ line duality. The result 
is a "line-curve" which is seen as the envelope of it's tangents. Usually 
this gives an unclear image and is at the heart of the "over-plotting" 
problem; a barrier in the effective use of Parallel Coordinates. This 
problem is overcome by a transformation which provides directly the 
"point-curve" representation of a curve. Earlier this was applied to 
conies and their generalizations. Here the representation, also called 
dual, is extended to all planar algebraic curves. Specifically, it is shown 
that the dual of an algebraic curve of degree n is an algebraic of degree 
at most n{n—l) in the absence of singular points. The result that conies 
map into conies follows as an easy special case. An algorithm, based on 
algebraic geometry using resultants and homogeneous polynomials, is 
obtained which constructs the dual image of the curve. This approach 
has potential generalizations to multi-dimensional algebraic surfaces 
and their approximation. The "trade-off" price then for obtaining 
planar representation of multidimensional algebraic curves and hyper- 
surfaces is the higher degree of the image's boundary which is also an 
algebraic curve in ||-coords. 
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1 Parallel Coordinates 



Over the years a methodology has been developed which enables the visu- 
alization and recognition of multidimensional objects without loss of in- 
formation. It provides insight into multivariate (equivalently multidimen- 
sional) problems and lead to several applications. The approach of Parallel 
Coordinates (abbr. ||-coords) p] is in the spirit of Descartes, based on a 
coordinate system but differing in an important way as shown in Fig^ On 
the Euclidean plane (more precisely on the projective plane P^) with 
j;y-Cartesian coordinates, n copies of a real line, labelled Xi, X2, . . . , X„, 
are placed equidistant and perpendicular to the x-axis with Xi and y be- 
ing coincident. These lines, which have the same orientation as the y-axis, 
are the axes of the Parallel Coordinate system for the n-dimensional Eu- 
clidean space M". A point C = (ci, C2, . . . , c„) € M" is represented by the 
polygonal line C having vertices at the values Cj on the Xi-axes. In this 
way a one-to-one correspondence is established between a points in and 
planar polygonal lines with vertices on the parallel axes. The polygonal line 
C contains the complete lines and not just the segments between adjacent 
axes. 

The restriction to provides that not only is a point represented by 



Figure 1: A point C = (ci, C2, . . . , c„) G M" is represented by the polygonal 
line C (consist of n — 1 sections) with vertices at the Cj values of the Xi axis 
for i = 1, 2, . . . , n. 
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a line, but that a line is represented by a point. The points on a line are 
represented by a collection of lines intersect at a single point as can be seen 
in Fig. [2 i.e. a "pencil" of lines in the language of Projective Geometry. 
A fundamental point ^ line duality is induced which is the cornerstone 
of deeper results in ||-coords. The multidimensional generalizations for the 
representation of linear p-flats in M" (i.e. planes of dimension 1 < p < n — 1), 
in terms of indexed points have been obtained. The proper setting for du- 
alities is the projective rather than the Euclidean plane M?. A review of 
the mathematical foundations is available in IHj. 




Figure 2: Fundamental point <-> line duality. 

For non-linear, especially non-convex, objects the representation is nat- 
urally more complex. In a point-curve, a curve considered as collection 
of points, is transformed into a line-curve; a curve prescribed by it's tangent 
lines as in Fig. |31 The line-curve's envelope, a point-curve, is the curve's 
image in ||-coords. In many cases this yields an image that is difficult to dis- 
cern. This point requires elaboration in order to motivate and understand 
some of the development presented here. As posed, the construction of a 
curve's image involves the sequence of operations: 

point — curve — > line — curve i — > 

point — curve {as the envelope of line — curve). 

Unlike the example shown in Fig. |3J where the line-curve's image is clear, 
the plethora of overlapping lines obscures parts of the resulting curve-line. 
This is a manifestation of what is sometimes called "over-plotting" problem 
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Figure 3: Point-curve mapped into line-curve (envelope of lines). 



in 1 1 -coords; the abundance of "ink" in the picture covers many underlying 
patterns. Simply our eyes are not capable of "extracting" the envelope of 
lots of overlapping lines. There are also computational difficulties involved 
in the direct computation of the image curve's envelope. In a way this is 
the "heart" of the "over-plotting" problem considered as a barrier in the 
effective use of || -coords. This problem can be overcome by skipping the 
intermediate step and go to an equivalent 



transformation which provides directly a clear image, as shown in Fig. [51 

This point-to-point mapping (Inselberg Transformation [3] ) preserves 
the curve's continuity properties. The idea is to use the — > part of the 
duality and map the tangents of the curve into points as illustrated in Fig. 
El Therefore, a curve is represented by a point-curve (the "dual curve") in 
the ||-coords plane. For a point-curve, 7, defined implicitly by 



the X, y coordinates of the point-curve image (i.e. dual) 7 are given by 



point — curve 



point — curve 



7 : f{xi,X2) = 0. 



(1) 



X 



df/dx2 



< 



{df/dxi+df/ax2) ' 



(2) 



y = 



{xidf /dxi+X2df /dx2) 
{df/dxi+df/dx2) 
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Figure 4: In line-curve (on the left) the tangents cover the image shown on 
the right as a point-curve. 



It was shown by B. Dimsdale ^ and generalized in [51, using this trans- 
formation that conies are mapped into conies in 6 different ways (Ij. Here we 
develop the extension of the dual image for the family of general algebraic 
curves. This family of curves is highly significant in many implementations 
and applications, since these curves can easily and uniquely be reconstructed 
from a finite collection of their points (for instance using simple interpolation 
methods (Tj, ^Sl)- Approximations of such curves can simply be obtained 
using similar methods. 

As will be seen, the dual of an algebraic curve in general has degree 
higher than the original curve. There are some fringe-benefits, illustrated 
later, where "special-points" such as self-intersections or inflection-points 
which are conveniently transformed. But we do not want to get ahead of 
ourselves. The key reason for this effort is to pave the way for the represen- 
tation of algebraic curves and more general hyper-surfaces, as well as their 
approximations, in terms of planar regions without losing information. To 
pursue this goal then we need firstly to study the image of curves starting 
with general algebraic curves. 
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Figure 5: Point-curve mapped into point-curve. 

2 Transforms of Algebraic Curves 
2.1 The Idea Leading to the Algorithm 

In order to represent non-linear relations in || -coords it is essential to extend 
the representational results first to algebraic curves; those described by ei- 
ther, explicitly or implicitly by irreducible polynomials of arbitrary degree. 
The direct application of eq. turns out to be difficult even for degree 
3. There is a splendid way to solve this problem using some ideas and tools 
from Algebraic Geometry. These involve properties of homogeneous poly- 
nomials and Resultant which are explained during the development of the 
method. For extensive treatments the reader is referred to [S], jUj and 

m- 

Starting with an algebraic curve 7 defined by an irreducible polynomial 
f{xi,X2) = in R2 its image in || -coords is sought. A number of prepara- 
tory steps smooth the way for the easier application of the transformations 
in eq. ((21) • First the curve 7 C is "raised" to a surface embedded in the 
projective space and only then the corresponding mapping into ||-coords 
is applied. 

Step 1 : There exists a one-to-one correspondence (preserving the reducibil- 
ity of polynomials) between any polynomial f{xi,X2) = of degree n 
and a homogenous polynomial in the projective plane P^, which is 



obtained by using homogeneous coordinates. Specifically, 

• replace 2;^ by ^ for /c = 1, 2, 

• multiply the whole polynomial by X3" and simplify. 

These multiplications are, of course, allowable for 7^ 0. The result 
is a homogeneous polynomial (which is also irreducible if / = is) 
with each term having degree n. To wit, 

n 

f{xi,X2)= ^ aijx\xl 

i+j=0 

n 

F{xi,X2,X3) = ^ aijx\x^2^l^~'~^ , (3) 

i+j=0 

which describes a surface in where the original polynomial curve is 
embedded as: 

/(X1,X2) =F(X1,X2,1). (4) 



Step 2 : The Gradient of F is found and denoted by : 

^r./ ^ ^dF dF dF^ . . M 

VF(xi,X2,X3) = (— , — , — ) = [r],LW)- 
oxi 0x2 0x3 

The three derivatives provide the direction numbers of the normal 
to the tangent plane at the point (xi,X2,X3) . It is a fundamental 
property of homogeneous polynomials that 

dF dF dF ^ 

xi^ hX2^ hX3-— = nF (5) 

oxi 0x2 0x3 

where n is the degree of F. In our case F = and hence the equation 
of the tangent planes is 

rjxi + S,X2 + ipX3 = 0. (6) 

Since the tangent plane at any point of the surface goes through the 
origin it is clear that the surface F is a cone with apex at the origin 
as shown in Fig. 121 



7 




Figure 6: Cone F{xi, X2, xs) 
f{xi,X2) = in p2. 



generated by embedding the curve 



Step 3 : Substituting for X3 from eq. © in 



F{xi,X2, 



,r]xi + ^X2 



)) = 0. 



(7) 



provides the intersection of the tangent plane with the cone which is 
a whole line as shown in Fig. [7| Each one of these lines, of course, 
goes through the origin and therefore it can be described by any one 
other of its points and in particular by (xi,X2,l). Simplifying the 
homogeneous coordinates by, 



{X1,X2, 



, r]Xi + Cx2 , 



{■ipxi, 11^X2, -{rjxi + ^3:2)) 



results in 



F(V'Xi, ■0X2, -{rjxi + ^X2)) = 0. 



(8) 



By the way, this is also a homogeneous polynomial in the five variables 
appearing in its argument. It is helpful to understand the underlying 
geometry. Eq. ^ can be considered as specifying a family of lines 
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Figure 7: The intersection of the cone F{xi,X2,X3) = with any of it's 
tangent planes is a whole line. Along such a line the direction numbers 
(??, ^, ip) are constant and unique. 



through the origin (0, 0, 0) and through each point {xi,X2, 1) along the 
curve given by eq. Q). Alternatively, the cone can be considered as 
being generated by a line, the generating line - see Fig. |H1 pivoted at 
the origin and moving continuously along each point (xi, X2, 1) of the 
curve described by eq. Q. On each one of these lines the direction 
numbers (ry, ^, ■0) are unique and constant as shown in Fig. El i.e. there 
is a one-to-one correspondence: 

Eq. 0, and hence it's rephrasing eq. (jHI, contains two equivalent 
descriptions of the cone, i.e. in terms of the coordinates (xi, X2, 1) and 
also the direction numbers (7^,^,-0). Hence one can be eliminated and 
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as it turns out it is best to eliminate the Xi , i = 1,2 something which 
is very conveniently done by means of the Resultant. The resultant 
R{F, G) of two homogenous polynomials 



F{xi,X2) 
G{xi,X2) 



and 



Er=o^i2;i^X2'"-^ =0 



is the polynomial obtained from the determinant of their coefficients 
matrix: 



R{F, G) 



( ao ai a2 



Det 



ao 



6o 

Vo 



hi 
b2 



bml 



where the empty spaces are filled by zeros. There are m lines of and 
n lines of hj. When F and G are both irreducible so is their resultant. 
Due to the homogeneity of F and G, 

R(F, G)=04^F = 0, G = 0. 

Let us rewrite eq. (|H| as 



i=0 

Appealing again to the homogeneity of F to obtain the relation 

dF dF 
xi- Vx2-^ = nF 

OXi OX2 



(9) 



which has already been mentioned earlier. Therefore 

F = 



axi 0x2 



and from the property of the resultant of homogeneous polynomials 
mentioned above 

dxi ' dx2 dxi ' dx2 
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Altogether then 

dF OF 

has finite degree since F and R are polynomials. It was shown that 

R = 0^ F = 



so not only their zero sets agree, which ensures that they are equiva- 
lent polynomials, but they also agree at an infinite number of points. 
Hence, R = provides the description of the cone F = in terms of 
the direction numbers (ry, ^, ip) and is also a homogeneous polynomial. 
The degree of ^ is n — 1. This and the structure of the resultant 
with the four triangular zero portions results in R being a polynomial 
of degree at most n[n — 1). 




Figure 8: This shows the generating line of the cone (which is found from the 
intersection of the tangent planes with the cone) together with the gradient 
vector. 
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Step 4 : Since we are interested in the solutions of = the multipher in 
powers of ip of R, if there is one, can be safely neglected due to the 
homogeneity. 



Step 5 : Only in this, the final step, the transformation of the curve from 
the xiX2-plane to the xy-plane with parallel coordinates is performed. 
It is done using equations Q rewritten, in view of eq. © as 

- « • (11) 



It is important to notice that both this and eq. (|1U|) involve only the 
derivatives r], ^, ip. Let c = r/ + ^, so that ^ = cx, rj = c — cx = c(l — x) 
and if) = —cy. Substitution provides the transform of the original 
curve /(xi, X2) = : 

R{V, V') = -R(c(l - x),cx, -cy) = 

c"'R{{l-x),x,-y) = 

R{{l-x),x,-y)=0, 

when c 7^ 0. In the previous step it was pointed out that this polyno- 
mial has degree at most n(n — 1), the actual degree obtained depends 
on the presence of singular points in the original algebraic curve / = 0. 



2.2 Algorithm 

T^ere the process involved is presented compactly as an algorithm i3_ whose 
input is an algebraic curve 7 : f{xi, X2) = and the output is the polynomial 
which describes 7, the curve's image in |[-coords. To emphasize, the algo- 
rithm applies to implicit or explicit polynomials of any degree and curves 
with or without singular points. The formal description of the algorithm is 
followed by examples which clarify the various stages and their nuances. 

For a given irreducible polynomial equation (otherwise apply the algo- 
rithm for each of it's component separately) f{xi,X2) = 0. 

1. Convert to homogeneous coordinates to obtain the transformation of 
/, a homogenous polynomial F{xi, X2, x^) = 0. 

2. Substitute 

xi ipxi, X2 — > ipX2 and 
X3 -{11x1 +CX2). 
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3. Find the resultant of the two derivatives F^^ and F^^ ■ 

4. Cancel the multiplier in a power of ip of the resultant R and denote 
the result by R'. 

5. The output is obtained by the substitution 

r) = 1 — X, i = X, ip = —y in R'. 

3 Examples of Algebraic Curves and their Trans- 
forms 

3.1 Conic Transforms 

The algorithm is illustrated with some examples starting with the conies. 

/ Ai A4 ArA / xi\ 

f{xi,X2) = {Xi X2 1) ^4 A2 ^6 = 

\A, Ae AsJ \ l J 

Aixj + 2^4x1x2 + 2^5x1 + A2XI + 2^63:^2 + ^3 = 0. 
Applying the algorithm in the sequence given above: 

Step 1 : The homogeneous polynomial obtained from / by using homoge- 
neous coordinates is 

F(xi,X2,X3) = 

Aixj + 2A4X1X2 + 2A5X1X3+ 

A2XI + 2^6X2X3 + ^3X1 = 0. 

Step 2 : Substituting 

Xi^lpXi, X2^iIjX2, ^3 ^ -(77X1 +^X2) 

yields 

F{'4>xi,'>pX2, -(r/xi + ^X2)) = 

(-2A5V'r? + AiV'^ + A^T]^) + xi^ 
2(-A5'0C + A^^ri + A^i^^ - A6i^ri)xiX2+ 
{A2^^ - 2Ae^^ + Ase)x2'' = 

CiXi^ + 2C2X1X2 + 2C3X2^ = 0, 

where q = Ci{r],(,,ip), for i = 1,2,3. 
Step 3 : Calculate the two derivatives of F and their resultant: 
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§^ = 2ciXi + 2C2X2, 
^ = 2C2X1 + 2C3X2, 

^dxi^dx2' \^2C2 2C3 y 

#^[(^3^2 - Al) r?2 + (A1A3 - Al) e+ 
{A1A2 - Al) + 2(^A6 - A^Ai)r]i+ 
2(A^4 - AiA^)iil^ + 2(^4^6 - 2^^2)#- 



Step 4 : Retain the resultants component which is multiplied by a power of 
ip and let 

Step 5 : The dual of /(xi, 2:2) = is then given in matrix form by 
R'{1- x,x,-y) = 

(oi 04 05 
04 a2 ae 
05 ae 03 

where the individual a,- are : 



02 = A1A2 - A42, 

03 = ^2^3-^6^, 

04 = ^6(^1 + Ai) - A^{A2 + A4), 
as = + - yl3(^2 + ^4), 

06 = A2A - ^4^6- 

The result is illustrated for the ellipse in fig Inland can also be contrasted to 
the earlier ways for obtaining the transformation 

eonie ^ eonie. 















-5 + ^6)' 



3.2 Algebraic Curves of Degree Higher than Two 

AAext the algorithm is applied to the algebraic curve of 3rd degree 

/(Xl, X2) = x\ - xl - x\ + X2 - I = ^■ 
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Figure 9: An ellipse is mapped into hyperbola. 




Figure 10: The cubic curve xf — — x| + X2 — 1 = (on the left) is mapped 
into a six degree curve (on the right). Notice also an instance of the duality 
inflection-point <-> cusp [4]. 



Step 1 : The homogeneous polynomial obtained from / by using homoge- 
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Figure 11: The self intersection point in the cure x\ + x^ — 8x1X2 = (on 
the left), disappears in the image curve 

270yx - 1987/x2 + AQyx^ + 108x^ - 36x^ - Blx^ = (on the right) . 
neous coordinates is 

F(xi, X2, X3) = X^ - xfxs - X2X3 + X2X| - x| = 0. 

Step 2 : Substituting 

xi^ipxi, X2 ^ ipX2, X3 ^ -(ryxi + ^X2) 
yields 

F{^pXl,^pX2,-ir]Xl +^X2)) = 
{xjj^ + ^'^r] — r]^)xi^+ 
(-3^7/2 + ^r]^ + ^2^)xi2x2 + 
{2ip^r] + tp'^r] - 3^^r/)xiX2^+ 

(-c^ + ve^ + v^o^2^ = 

CiXi^ + C2Xi^X2 + 03X1X2^ + C4X2^ = 0, 

where the q = Cj(r/, ^p), for i = 1, 3. 
Step 3 : Calculate the resultant of the two derivatives of F: 

§^ = SciXi'^ + 2C2X1X2 + C3X2^, 
W- = C2Xl^ + 2C3X1X2 + 3C4X2^, 
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/3C1 2C2 C3 \ 

^ — r)^+ 

dx\ ' dx2 



3ci 2c2 C3 

C2 2C3 3C4 



\ C2 2C3 3C4 / 

where i?' is polynomial of degree 6. 
Step 4 : Discarding the resultant's factor {—ip^) i.e. let 

Step 5 : Finally, substituting 

rj = 1 — X , ^ = X , ■0 = ~V 
results in the 6th-degree curve : 
R'{1 -x,x,-y) = 

-23 + 292y22;2 _ ^22x'^ + 326yx^ - U6y'^x 

+610x3 + 237/2 - 27yS2 + - 27y'' 

-22yx^ - 244y2x3 - - 420yx2 

-126y3x + 232yx + 214x^ - 499x'' + DOy^x^ 

+54y3 + 156x - 50y - 31x^ + 71^2^4 

-lAy^x^. 

The source and image curve are shown in Fig. IIUI In this case n = 3 and 
the maximum degree n(n — 1) = 6 is attained. Notice that this covers as 
special cases : 

- the point line duality where points (with n = 0) are mapped into lines 
(with n = 1) and vice- versa, and 

- the conies with n = 2. 

The next two examples show that the image curve may have degree less 
than n(n — 1). Apparently the maximal degree is attained by the image 
curve in the absence of singularities in the function or it's derivatives. The 
full conditions relating the image's degree to less then the maximal are not 
known at this stage. Two such examples are shown in the subsequent figures, 
Fig. HUand Fig. [H 

4 Conclusions 

Consider the class S of hyper-surfaces in which are the envelopes of their 
tangent hyper-planes. As mentioned earlier jH] hyper-planes can be repre- 
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Figure 12: xfx2 -1 = i — > + 27x^ - 27x2 = 0. 



sented in ||-coords by n — 1 indexed points. Hence for a hyper-surface a € 5 
each point P € o" maps into n — 1 indexed planar points. From this it fol- 
lows that the hyper-surface a can be mapped (i.e. represented) by n — 1 
indexed planar regions ai composed of these points. Restricted classes of 
hyper-surfaces have been represented in this way and it turns out that the 
(Ti reveal non-trivial properties of the corresponding hyper-surface a. It has 
already been proved that for any dimension Quadrics (algebraic surfaces 
of degree 2) are mapped into planar regions whose boundaries are conies 
[1]. This also includes non-convex surfaces like the "saddle". There are 
strong evidence supporting the conjecture that algebraic surfaces in general 
map into planar regions bounded by algebraic curves. This, besides their 
significance on their own right, is one of the reasons for studying the rep- 
resentation of algebraic curves. In [2] families of approximate planes and 
flats are beautifully and usefully represented in ||-coords. Our results cast 
the foundations not only for the representation of hyper-surfaces in the class 
(5), but also their approximations in terms of planar curved regions. 
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